Epitaxial TbMnOa thin films on SrTiOs substrates: A structural study 



C.J.M. Daumont^, D. Mannix^, Sriram Venkatesan^, G. Catalan^, 
D. Rubi\ B.J. Kooi\ J.Th.M. De Hosson\ B. Nohedai 
^ Zernike Institute for Advanced Materials, University of Groningen, 9747 AG Groningen, The Netherlands 
'^Institut Neel, CNRS-UJF, 25 Avenue des Martyrs, 380^2 Grenoble Cedex 9, France and 
Earth Science Department, Gambridge University, Gambridge GB2 3EQ, United Kingdom 

TbMnOa films have been grown under compressive strain on (OOl)-oriented SrTiOa crystals. They 
have an orthorhombic structure and display the (001) orientation. With increasing thickness, the 
' structure evolves from a more symmetric (tetragonal) to a less symmetric (bulk-like orthorhombic) 

I structure, while keeping constant the in-plane compression thereby leaving the out-of-plane lattice 

, spacing unchanged. The domain microstructure of the films is also revealed, showing an increasing 

' number of orthorhombic domains as the thickness is decreased: we directly observe ferroelastic do- 

^ I mains as narrow as 4nm. The high density of domain walls may explain the induced ferromagnetism 

observed in the films, while both the decreased anisotropy and the small size of the domains could 
account for the absence of a ferroelectric spin spiral phase. 

PACS numbers: 68.55.-a 



O . I. INTRODUCTION 

Multiferroics are a class of materials exhibiting coexistence of two or more ferroic orders [H, i, HQ: (Anti- 
)ferromagnetism, (anti-)ferroelectricity, ferroelasticity and/or ferrotoroidicity. Multiferroic materials have been in- 
creasingly studied in recent years due to the possibility of substantial coupling between the ferroic properties within 
a single phase. The most interesting case for applications occurs when the magnetic and electric degrees of freedom 
are coupled, giving rise to a large magnetoelectric (ME) effect. This can eventually lead to a range of novel devices 
that use the control of the spontaneous magnetization (polarization) of a material with an electric (magnetic) field. 
Unfortunately, multiferroics are very rare and efforts have been made to discover and synthesize new multiferroic and 
magnetoelectric materials, as well as to understand the underlying mechanisms. 
Q ' Among the possible mechanisms producing multiferroicity and large magnetoelectric coupling in single phase ma- 
O I terials, the induction of ferroelectricity in spiral spin systems has been intensively studied [y, |7|. The orthorhombic 
perovskite TbMnOa is the best known example of such a system[3, H, Q and has also been shown to present giant 
' ME effect^. Among the rare earth manganites, TbMnOs presents an intermediate Mn-O-Mn bond angle and lies, 
^ in the temperature versus ionic radius phase diagram, in between two different magnetic phases, displaying a com- 
plex magnetic behavior|9j. TbMnOs shows antiferromagnetic ordering below 40K, where the Mn spins align 

in an incommensurate sinusoidal structure. By decreasing the temperature further, the propagation vector of the 
sinusoidal spin structure decreases until it locks at T/ocfc~ 28K, where the magnetic structure changes into an spiral 
. structure [sl-M [ol. [l0| that propagates along the b-axis. Due to the Dzyaloshinsky-Moriya(DM) interaction[ll|, a 
T-H ' spontaneous polarization along the c-axis and a strong ME effect arises below T/ocfe- 

For practical devices, multiferroics are often preferred in thin film form. Moreover, the strain induced by the 
mismatch between the film and the substrate lattice parameters can lead to improving the materials properties with 
respect to the bulk. Interestingly, although TbMnOsis by now a very well known material, only a few reports on 
this perovskite in thin film form are available [Tsl uB^ . Since the magnetic and ferroelectric phases are determined 
by the crystal structure, a thorough characterization of the structure of the films is needed. This, however, is not 
straight forward due to the small thickness of the films and grazing incidence diffraction, using synchrotron sources, is 
' required. In this work, we present the growth and structural characterization of TbMnOs perovskite thin films with 
thicknesses lower than lOOnm on single crystals of (001)-SrTiO3. (001)-SrTiO3 is preferred in studies of epitaxy in 
perovskites because it can be obtained with atomically flat surfaces, favoring high quality growth. The films are shown 
to be very flat, c-oriented and epitaxially strained. Their unit cell is orthorhombic with an orthorhombic distortion 
that can be tuned with the film thickness up to a thickness of about 50nm. 

Although ferroelectric and magnetoelectric measurements are still in progress, the magnetic properties (reported 
elsewhere[14]) suggest that these films are remarkably different from the bulk material: No evidence of the spiral spin 
structure that gives rise to multiferroic behavior was found in the magnetic data. On the other hand, ferromagnetic 
interactions are induced in the films, which renders them appealing for applications such as spin valves, where 
ferromagnetic insulators are required. This is in agreement with a very recent work on the same system[l5] and 
reports in other orthorhombic manganite films[16]. We show here that the structure of the films and microstructure 
are likely to explain the difference between the films and bulk behavior. 
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FIG. 1: (Color Online): (a)Intensity of the specular spot of the RHEED pattern as a function of deposition time. The inset 
is a blow-up of the first RHEED oscillation showing the laser pulses, (b) Typical RHEED pattern obtained after growth and 
annealing. The picture is taken after cooling down to room temperature and under vacuum, (c) and (d) AFM images of an 
2nm and 8nm thin film, respectively, grown at p02= 0.25 mbar. The vertical and horizontal scales are shown at the bottom of 
the pictures 

II. EXPERIMENTAL 

(OOl)-oriented TbMnOs (TMO) thin films were deposited on (001)-SrTiO3 (STO) substrates by Pulsed Laser 
Deposition (PLD) assisted by reflective high energy electron diffraction (RHEED), using a KrF excimer laser with 
A = 248nrn. A stoichiometric target of TbMnOa was prepared by means of a standard solid-state reaction. Prior 
to deposition, the substrates were chemically treated and fired in O2 in order to obtain Ti02 single-terminated flat 
terraces (it). The deposition of the films reported here was performed at 750°C in oxygen pressures of 0.25 mbar to 
0.9 mbar. The laser fluencc was 2 J/cm^, with a repetition rate of IHz, and a distance between substrate and target 
of 55mm. After deposition, the films were slowly cooled down (-3°C/min) to room temperature under an oxygen 
pressure of 100 mbar. 

The highly distorted perovskite TMO crystallizes in a orthorhombic structure (space group: Pbnm) with lattice 
constants of ao=5.2931 A, bo=5.8384 A and Co=7.4025 A [ISi]- However, lattice parameters reported so far for the 
bulk/single crystal system, and in particular the c-lattice parameter, vary significantly and this has been attributed 
to different growth conditions 0, 0, [lO]. The STO substrate is a cubic perovskite with lattice constant of a=3.905 A. 
During growth, the intensities of the RHEED patterns change as a result of the relative surface coverage, roughening 
and disorder of the growing layers [2lj. The intensity of the specular spot of the RHEED pattern was recorded from 
the beginning of the growth (see Figure 1(a)). With the first laser pulses, a large drop of the RHEED intensity takes 
place, indicating an important roughening of the surface with respect to the STO substrate. The recovery of intensity 
observed after approximately 35 laser pulses corresponds to the growth of the first complete layer of TMO. 

The crystallinity and structure of the films were studied by standard x-ray diffraction (XRD) with a Panalytical 
X'Pert MRD diffractometer, while their surface morphology was probed by means of a Nanoscope Ilia atomic force 
microscope (AFM). X-ray photoelectron spectroscopy (XPS) experiments were performed in a SSX-100 spectrometer 
from Surface Science Instruments with a monochromatic Al-K x-ray source (hz^=1486.6eV). XPS in films of different 
thickness and grown at different oxygen pressures are consistent with a manganese valence of +3 and no evidence of 
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Mn4+ or Mn2+ was found [T3|- In order to probe the in-plane lattice, grazing incidence x-ray diffraction (GIXD) 
was performed with synchrotron radiation, both at the Wl beamline in HASYLAB-DESY (Hamburg) and the XMaS 
beamline in ESRF (Grenoble). TEM plan- view specimens were prepared by conventional procedure involving cutting, 
grinding, polishing, dimpling and ion milling. A precision ion polishing system (Gatan model 491) with 4 Kv Ar-|- 
beams was used, in which both guns make an angle of 8° with the bottom side of the substrate. On the other side, 
the film is covered with a piece of glass to avoid any contamination by redeposition from the sputtered backside. The 
observations have been performed with a JEOL 201 OF electron microscope at an accelerating voltage of 200 kV. 



III. RESULTS 



About eight RHEED intensity oscillations can be seen at the beginning of the growth corresponding to an initial 
layer-by-layer (2D-like) growth. However, those oscillations are superimposed on a decreasing intensity background, 
indicating that an overall roughening occurs as the growth proceeds. After the oscillations have faded away and the 
intensity of the specular spot has decreased to a certain intensity, the later remains constant during the rest of the 
deposition (outside of the graph range) and the pattern changes from a purely stripy pattern to the mixed stripe-spot 
pattern shown in Figure 1(b). This type of behavior has been reported recently for other oxides as a 'pseudo-2D islands 
growth' [13, [13 ■ This change from 2D character to 'pseudo'-3D character of the growth has also been identified using 
AFM. Indeed, Figure 1(c) shows the morphology of a 2nm film, for which the deposition was stopped shortly before 
the 2D-to-3D transition. It can be seen that the film is atomically flat, showing the step- and- terrace morphology of 
the substrate. However, when the deposition is stopped at the initial stage of the transition from 2D to 3D growth, 
particles with a diameter of around 30 nm start appearing, as shown in Figure 1(d). We found that the change from 
2D to 3D growth occurs at a thickness in between 2 and 5 nm. 
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FIG. 2: (Color Online): AFM images of 17nm (a), 34nm (b) and 67nm (c) thick TMO thin films grown with a p02= 0.25 
mbar. The scans area is l^mxl/xm. The corresponding diffraction patterns can be seen in (d), (e) and (f), respectively. The 
experimental data are are shown as thick black lines. The fit to the data is shown in each case as a thin blue line. 

Figures 2(a) to 2(c) show AFM images of three TMO films with different thicknesses grown under the same 
conditions. In the 17nm film (Figure 2(a)), although it is already in the 3D growth regime and the grains dominate 
the film morphology, the steps from the substrate are still visible. The steps are not clearly visible for a 34nm film 
and a 67 nm (Figures 2(b), (c)). Figures 2(d)-(f) show typical 2d-uj x-ray diffractograms around the (OOl)c Bragg 
reflection of STO (the most intense ones in the pattern), for the same three films. The film peaks are those at the 
right hand-side of the STO refiections. The blue lines are the fits of the crystal truncation rod using a kinematical 
model. These calculations take into account, not only the complex refractive index of the substrate and the film, but 
also the angular dependent atomic scattering factors. The high quality and fiatncss of the interfaces is evidenced by 
the Lane fringes around the film peak for thicknesses up to 67nm. 

Figure 3 shows a broader 28-uj scan including both the (OOl)c and the (002)^ STO Bragg reflections. No secondary 
or impurity phases could be observed. A phi-scan around the (024)pc (ie (228)o ) of TMO (see inset) shows the peaks 
separated 90 degrees from each other, confirming the four-fold symmetry of the film and the coherent growth. From 
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the position of the film peaks of figure 3, an out-of-plane lattice spacing of 3.737 A was found. X-ray diffraction, thus, 
shows that the films are single phase and (OOl)-oriented (with the c-axis perpendicular to the surface). Using a value 
of Poisson ratio of 0.35, typical in manganites[2l|. the out-of-plane lattice spacing is estimated to expand from 3.70 
A (c/2 of bulk TMO) to about 3.74 A, in good agreement with our measurements. 
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FIG. 3: 2e-uj for a 67nm thick film of TMO grown on STO. Inset; Phi-scan around the (228)o reflection. 

A larger out-of-plane lattice parameter than the bulk's one could also be explained by a film whose strain is relaxed 
by means of oxygen vacancies, as they are known to expand the unit cell[25l. |26| . To investigate this possibility, 
the oxygen pressure during deposition was increased up to 0.9 mbar (the maximum attained in our set-up). Within 
the range of oxygen pressure investigated, the out-of-plane lattice parameter indeed decreases linearly with increasing 
oxygen pressure during deposition, consistent with the number of oxygen vacancies decreasing for increasing pressures. 
The maximum pressure of 0.9 mbar seems to be large enough to produce stoichiometric films since the unit cell of the 
relaxed layers is as large as the bulk's one, as it will be discussed later. 
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FIG. 4: (Color Online): Reciprocal space maps around the (103)c Bragg reflection of the STO substrate for 17nm, 34nm and 
67nm films grown at 0.25 mbar. The vertical and horizontal lines are guides to the eyes. The abscissa (ordinate) represents 
the in-plane (out-of-plane) component of the scattering vector. Both are normalized by 2ko=47r/A. 



A detailed investigation of the film structure was performed using x-ray diffraction to map selected areas of the 
reciprocal space. Figures 4(a)-(c) show reciprocal space maps (RSMs) around the (103)c Bragg reflections of a 17nm, 
34nm and 67nm thick film, respectively, grown at 0.25mbar. The peaks corresponding to the films (substrates) are 
those at larger (smaller) Kpe^p (see horizontal lines in the figure). The in-plane and out-of plane components of the 
(103)c and (013)c TMO peaks do not change as a function of thickness and the in-plane lattice parameter remains 
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identical to that of STO for thicknesses up to 67nm, indicating that the films are coherent with the substrate along 
the pseudo-cubic [100] or [010] directions. As the thickness is increased, only a broadening of the rocking curve of 
the TMO film can be seen, indicating increase of mosaicity. However, when the RSM is taken around the (113)c, a 
different picture is observed, as shown in figure 5(a) for a 67nm thin film. Instead of a coherent film peak, two film 
peaks are found that could be indexed as orthorhombic (206)o and (026)o. The data therefore, fit a structural model 
in which the films are orthorhombically distorted, similar to the bulk material, but keep coherence with the cubic 
substrate along the [100] or [010] directions. 




FIG. 5: (Color Online): (a)RSM around (206)o and (026)o for a 67nni thick TMO film grown at 0.25 mbar. The arrow indicates 
the kpar of the STO substrate. (b)RSM around the (103)c STO Bragg reflection for a 80nm thick TMO film grown at 0.9 
mbar. The abscissa (ordinate) represents the in-plane (out-of-plane) component of the scattering vector, both normalized by 
2k<,=47r/A. 



The orthorhombic Aq, bg and lattice parameters can be extracted from the (116)o and (206)o reflections shown 
above. Figure 6 shows the evolution of the lattice parameters as a function of thickness for the films grown with 
p02= 0.9 mbar and p02= 0.25 mbar. Both series display similar trend, showing that the oxygen vacancies, which 
are likely to be present in the films grown at the lowest pressure, giving rise to larger Cq values, do not change the 
general behavior. The orthorhombic distortion is the smallest for the thinnest films and a 2nm thick TMO film shows 
a tetragonal symmetry as evidenced by GIXD (not shown here). With increasing thickness, Co remains unchanged, 
while Oo decreases and bo increases towards their bulk value. However, the bulk lattice parameters are not reached in 
a continuous way and Oo and bo saturate for a thickness of around 50 nm. Above about 7Gnm, a new spot appears in 
the diffraction maps, as shown in Figure 5(b). This second spot corresponds to a relaxed unit cell with bulk lattice 
parameters (dashed lines in Figure 6). The fact that, in the films grown at 0.9mbar, this spot appears with lattice 
parameters almost identical to the bulk ones, strongly indicates that these films are stoichiometric, as discussed above. 
We have shown that the films are orthorhombic but keep coherence with the cubic substrate along the [010]c or [100] c 
directions. Thus, the films orthorhombic lattice parameters fulfill + b^ = {lasroY for all thickness and only the 
in-plane angle 7 (see figure 9) changes with increasing thickness. Therefore, the strain state is constant and so is the 
c lattice parameter. 

In order to better understand the in-plane coherence of the films with the substrate and the domain structure of the 
film, GIXD experiments were performed using synchrotron radiation. Figure 7 (a)-(d) shows in-plane RSM's around 
selected reciprocal lattice points of the STO lattice for a 8 nm thin film grown at 0.25mbar. In all the maps, the 
central sharp spot corresponds to the STO Bragg peak and the other broader visible reflections around it are from 
the TMO film. As shown in (a), four broad peaks can be seen around the (110) substrate spot. They correspond to 
the different equivalent orientations of the orthorhombic unit cell (see sketch in figure 9). The presence of these four 
domains maintains in the film the four-fold symmetry of the cubic substrate. Moreover, as shown in figure 7(b) and 
(c), a peak broadening in the perpendicular direction (along the [010] for the (100) reflection and along [100] for the 
(010) reflection), the so-called rocking curve, is seen around the central cubic position. This corresponds to the part 
of the film that is coherent with the substrate (with lattice spacing 3.90 A). As seen in Figure 7(d), around the (220)c, 
only two of the four peaks can be observed, those corresponding to the (400)o reflections. The absence of the (040)o 
peaks is due to the relative intensities of those reflections. Indeed, the intensity ratio I(200)o/I(020)o is 1.4 in bulk, 
whereas a ratio of 77 is found for I(400)o/I(040)o. This makes the intensity of the (040)o too weak for us to detect it. 

GIXD experiments were also performed for a 55nm fllm grown at 0.25mbar. Figure 8 (a)-(d) shows different in- 
plane RSMs for this thicker fllm. From (a), only two film reflections are observed around the substrate (110) peak 
(which in this case is not visible due to the grazing incidence geometry and the larger thickness of the film). The 
peaks correspond to the (200)o (higher H and K values) and the (020)o (at lower H and K values). Analysis of these 
diffraction maps indicates that, along with the change of s.o and bo, a gradual rotation of the unit cell also occurs at 
larger thickness, so that the orthorhombic axes {&o and bo) align parallel to the [110] of STO (see figure 9c). This can 
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FIG. 6: Orthorhombic lattice parameters (flo, bo, Co/V^) as a function of thickness for films grown with p02= 0.9 mbar (filled 
symbols) and 0.25mbar (open symbols). In the 80nm film, two phases are found and one of them (symbols marked with a star) 
corresponds to the bulk orthorhombic unit cellfisl] . The lines are guides to the eye. 




FIG. 7: (Color online): RSMs of a 8nm TMO film grown at 0.25mbar around: (a) the (110)^; (b) the (010)^; (c) the (lOO)c 
and (d) the (220)c STO Bragg reflections. The axes are in reciprocal lattice units of the substrate (1 r.l.u.= 27r/3.905 A) 



be better seen in figures 8 (a) and (d), where the extremely broad peak is due to domains existing at aU orientations 
in between those sketched in figure 9(a) and (b) and those in (c). The coherence with the substrate in the thinner 
fihus is then given by [110]o parallel to [100]c (or to [010]c). By further increasing the thickness, the strain energy has 
increased such that the coherency along the [100] c cannot be kept and the films relax to their bulk structure, with 
the orthorhombic in-plane axis parallel to the [110] of STO. 

Figure 10 shows TEM plan view images and electron diffraction patterns recorded for two samples with thicknesses 
of 17nm and 67nm. In the plan view image of the 17nm TMO film, two differently oriented types of rectangular 
domains can be seen along with strain contrast for domains having the same orientation. This is consistent with the 
four types of domains seen via XRD. For the thicker film, two orientations can still be seen along with bigger domains, 
consistent with the coalescence of the grains seen by AFM. Moreover, strain contrast is present within the domains 
showing that the domains coalesce as the thickness increases. Screw dislocations can be seen on the thicker film. In 
addition of screw dislocations at the interface, another type of defects probably exist at the domain boundaries that 
take care of the in-plane rotation of the domains. The electron diffraction patterns confirm what has been seen in 
x-ray diffraction: for the thinner film, the diffraction peaks are the superposition of four TMO domains along with 
the contribution from the substrate; for the 67nm TMO thin film, the electron diffraction pattern differs from that of 



FIG. 8: (Color online): RSMs for a 55nm TMO film grown at 0.25mbar around (a) the (llO)c; (b) the (OlO)c; (c) the (lOO)c 
and (d) the (220)c STO Bragg reflections. The axes are in reciprocal lattice units of the substrate (1 r.l.u.= 27r/3.905 A) 
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FIG. 9: (a) and (b)Sketch of the four types of orthorhombic domains present in the thinner films and their coherence with 
the substrate. (c) Sketch of the orthorhombic domains present in the thicker films and their coherence with the substrate, d) 
Relationship between the orthorhombic and pseudo-cubic lattice parameters with the substrate lattice 



the very tliin films and shows a rotation gradient, as discussed above. 



IV. DISCUSSION 



For the 17nm film, the domains are found to have a small width of about 4nm (see Fig. 10a). This domain size 
is particularly small when compared with the typical thickness of magnetic domain walls, which are of the order of 
several nm[27','28'|. This means that in the thin films the fraction occupied by domain walls is comparable to that of the 
domains themselves, and therefore the properties of the domain walls are likely to affect the overall magnetic properties. 
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that is, the induced weak ferromagnetism, and perhaps even the absence of a lock-in transition [lj|. For example, 
while the orthorhombic domains are known to be antiferromagnetic in character, in agreement with the negative 
Curie- Weiss temperature measured in the films the magnetic interactions at the domain walls can give rise to a 
ferromagnetic component. Symmetry arguments show that magnetoelectric coupling can induce ferromagnetism in 
the domain walls of ferroelectric antiferromagnets[2^], and works by Fiebig and co-workers have also shown that the 
ferroelectric domain walls of multifcrroic hexagonal manganites can have a net magnetization at their center [sO] as 
well as enhanced magnetoelectric coupling [sij. 

The small domain size is likely to also affect the spin spiral. In bulk TbMnOs, the period of the spiral is of the 
order of 2nm, which is comparable to the domain size and so it may well be destroyed by the proximity of the domain 
boundaries. Moreover, even if the spiral spin order that gives rise to the lock-in and ferroelectric transition survived, 
the domain size may be inconmensurate with the spiral period, meaning that some spins will not be fully compensated, 
an effect that will be more noticeable for smaller domains. The presence of small domains would also prevent long 
range coherence, since the spiral propagation direction, which is along the orthorhombic b-axis, has to to flick from 
X to y directions at each domain wall. 

Finally, the reduction of in-plane anisotropy means that the difference between the a and b lattice parameters in 
the orthorhombic structure decreases. For the extreme case of the tetragonal films (2nm thick or less) there is no 
privileged direction, and hence there can be no spin spiral. But even for the thicker orthorhombic films, strain is likely 
to affect, or even destroy, the spin spiral, since the changes in in-plane anisotropy with respect to the bulk compound 
must also affect the Mn-O-Mn bond angle, which is known to be the critical parameter for the appearance of the 
lock- in spiral phase (3^ 




FIG. 10: Plan view TEM pictures of a 17nm (a) and 67nm (b)TMO film. Electron diffraction patterns of the same films are 
sliown in (c) and (d), respectively. The white circles shows show the presence of dislocations. 



V. CONCLUDING REMARKS 

We have successfully deposited epitaxial orthorhombic TbMnOs films on (001)-SrTiO3, with thicknesses ranging 
from 8nm to 80nm. The crystal structure of the thin strained films has been identified as a less distorted orthorhombic 
unit cell compared to the bulk one with &fiim>^buiki i>fiim<^buik and Cfiim>Cbuik- We found that the films grow 
with the c-axis out of plane and they orient themselves in the plane such that the (llO)o directions align with the 
(100}c directions of the substrate. The orthorhombic lattice parameters Qo and bo are constraint by epitaxy so that 
a,o^+bo^= (2a5To)^- This allows for four equivalent orientations and, therefore, four types of domains are observed, 
such that the in-plane diffraction patterns display four-fold symmetry. 
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With increasing thickness, the in-plane orthorhombic axes gradually change towards their bulk values, increasing the 
orthorhombicity of the films, but still keeping the partial coherence and epitaxy relation ao^+bo^= {2asTo)^- Due to 
this, the amount of in-plane compression (and hence out-of-plane elongation) remains unchanged up to thicknesses of 
about 70nm. Above this thickness, the lattice relaxes to the bulk one. In spite of the unchanged in-plane compression, 
the thicker films do show a distribution of unit cell rotations coupled with increasing orthorhombic distortion. 

It has been shown [l4[ that the physical properties of these films are very different from those of the bulk: The films 
show ferromagnetic interactions and spin-glass-like behavior below the Neel temperature of ~ AOK, which are absent 
in the bulk. This can be of interest due to the scarcity of ferromagnetic insulators. Similar results have recently been 
obtained for this and other orthorhombic perovskitesjisl IToj . Here, we propose that this is directly associated to the 
domain walls, showing evidence that these constitute a large part of the volume of the films. On the other hand, the 
magnetic anomaly observed in the bulk material at about 28K, related to the transition to the spiral spin structure 
and the onset of ferroelectricity, could not be observed in these films [l4|. This may be due to the decrease of in-plane 
anisotropy that we show takes place with decreasing thickness, or to the small domain size. Work is in progress to 
investigate these possibilities. 
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